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We make use of the formaUsm developed in Ref. []|, and calculate the chargino 

(N 

^ ' mediated baryogenesis in the Minimal Supersymmetric Standard Model. The formal- 

cn ■ 

. ism makes use of a gradient expansion of the Kadanoff-Baym equations for mixing 

•/^ \ fermions. For illustrative purposes, we first discuss the semiclassical transport equa- 
O 

' tions for mixing bosons in a space-time dependent Higgs background. To calculate 

r"| ' the baryon asymmetry, we solve a standard set of diffusion equations, according 

' , to which the chargino asymmetry is transported to the top sector, where it biases 

(— I \ sphaleron transitions. At the end we make a qualitative and quantitative comparison 



of our results with the existing work. We find that the production of the baryon 
asymmetry of the Universe by CP-violating currents in the chargino sector is strongly 
constrained by measurements of electric dipole moments. 

PACS numbers: 98.80.Cq, ll.30.Er, 11.30.Fs 



I. INTRODUCTION 

Electroweak baryogenesis Q] is an effective framework for explaining the baryon asym- 
metry of the universe (BAU). The most appealing feature of this mechanism lies in the fact 
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that the relevant physics will soon be explored by experiments, most notably by LHC at 
CERN and by the new generation of electric dipole measurements. 

It has been realized that the scenario of electroweak baryogenesis depends on extensions 
of the Standard Model (SM), since two mandatory conditions are not met in the SM. The 
first reason is that CP-violation in the SM is marginal, such that the observed magnitude of 
baryon asymmetry cannot be explained. Secondly, the electroweak phase transition in the 
SM is a crossover ^, |^, leading to a too weak departure from equilibrium to be viable for 
baryogenesis. 

The Minimal Supersymmetric Standard Model (MSSM) instead has all the necessary 
ingredients. CP violation is enhanced by adding phases to the parameters in the soft su- 
persymmetry breaking sector, which contribute to the chargino mass matrix. Furthermore, 
the additional bosonic degrees of freedom can lead to a strong first order phase transition 

These considerations indicate that the MSSM has the potential of explaining the observed 
BAU via electroweak baryogenesis. However, a formalism that determines the baryon asym- 
metry has to incorporate several features. Clearly the formalism has to reflect the quantum 
nature of the involved particles, for CP violation is a purely quantum effect. In addition, 
since the sphaleron processes are only operative in the unbroken phase, the CP-violating 
particle densities have to be transported away from the wall into the unbroken phase to 
lead to a net baryon density. A formalism that can handle both of these aspects is given 
by the Kadanoff-Baym equations, which are in turn derived from the out-of-equilibrium 
Schwinger-Dyson equations. 

Early approaches that aimed to determine CP- violating densities and have not attempted 
to derive transport equations from first principles have been based on the dispersion rela- 
tion of the quasi-particles 0, Isjs, ll^ llll deduced with the WKB method. For a recent 



resurrection of the method see ll 
In 



ecui 





IJ] it was suggested that an important contribution is given by mixing effects 
of the quasi particles in the wall rather than from the dispersion relations in the case of a 



gisp ersi' 

HQ 



nearly degenerate mass matrix. However, in the work [Ij, transport equations are not 
derived in a first principle approach either, but the current continuity equation is used to 
determine CP-violating contributions to the Green functions in a perturbative approach, 
which are subsequently inserted as sources into classical diffusion equations derived in 
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These classical diffusion equations neglect oscillations of the off-diagonal elements of the 
Green function that are important for a proper treatment of CP violation. 

Starting from the Kadanoff-Baym equations, the authors of j^, Q] have derived the 
CP-violating semiclassical force in kinetic transport equations, which appears in fermionic 



kinetic equation at second order in derivatives. Initially, this was done for the one fermion 
flavour case and then subsequently generalized to the diagonal part of the multifiavour 



case 11^ 

Recently this formalism was advanced to include mixing fermions Q|. The formalism 
provides an accurate description of the dynamics in the thick wall regime, which applies 
to particles, whose de Broglie wave length is much shorter than the thickness of the phase 
boundary (bubble wall), formally <^ k. 

lat two features of the transport equations are not 



One conclusion of the work [1| is t 
captured by the procedure used in [13, 



14| . Firstly, the densities that are off-diagonal in the 



mass eigenbasis of the system will perform oscillations analogously to neutrino oscillations. 
This effect suppresses the transport of the CP-violating sources, especially if the mass spec- 
trum in the chargino sector is far from degeneracy. Secondly, while Refs. used a 
phenomenological prescription (Pick's law) to introduce the CP-violating sources into the 
diffusion transport equations, no such prescription is required in our formalism. The sources 
enter the diffusion transport equations with an unambiguously defined amplitude. 

A first goal of this publication is to study the simpler bosonic case, and thus to rectify 
the conclusions of ^. As a second and principal goal, we consider the chargino mediated 
baryogenesis in the MSSM, in order to study the effects of flavour oscillations and source 
amplitude ambiguity on the baryon asymmetry within the framework of the reduced set 
of diffusion equations for charginos and quarks used in We also make a 

comparison of baryogenesis from the semiclassical force mechanism. The principal difference 
with respect to the previous work, is that our treatment is basis independent, while the 
calculations presented in Ref. |lO, ll2l ll^ were performed in the mass eigenbasis. 

This article is organized as follows. In section |n] we derive transport equations for mixing 
bosons. This is done mainly to clarify the conclusions from [1] that are present in the bosonic 
case, too. In the subsequent section we discuss, how the introduction of phenomenological 
damping terms can lead to additional unphysical CP- violating sources. The sections HVl and 
Instate the fermionic transport equations derived in Q and the system of diffusion equations 
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that is used to determine the baryon asymmetry. Numerical results are presented in section 
IVH and we conclude in section lyTTl 



II. TRANSPORT EQUATIONS FOR MIXING BOSONS 

In this section we will derive transport equations for mixing bosons from the Kadanoff- 
Baym equations and the resulting CP violating particle densities. This is a simpler analogon 
to the derivation for the fermionic case given in |lj]. In the fermionic case, the spinor structure 
complicates the decoupling of the system of equations, but the bosonic case given here will 

n 

already support the main conclusions given in y^] without the technical issues coming from 
the spinor structure. 



A. Kadanoff-Baym Equations and the Approximation Scheme 

Starting point are the coupled Kadanoff Baym equations 

e-Ojp _ M^}{A<^>} - e~'^{Uh}{A<'>} - e-^^{n<'>}{AJ = C, (1) 
C = ie-0({n>}{A<}-{n<}{A>}), (2) 

where A denotes the Green function and 11 the self-energy of the bosons. Both quantities 
are N x N matrices in flavour space and depend on the average coordinate and the 
momentum variable fc^. The superscripts <, > and the subscript h denote the additional 
2x2 matrix structure as usual in the Kadanoff-Baym formalism 




A< = A+-, A> = A-+, A* = A++, A* = A"-, 

A, = A*-|(A< + A>). (3) 

The diamond operator coming from the transformation into Wigner space is defined by 

0{a}{b} ^ \ [{dx,a)d,h - {d,^a)dxh) . (4) 
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The mass squared matrix is space-time dependent and hermitian. During the elec- 
troweak phase transition, the bosonic particles relevant for baryogenesis are the squarks 
whose mass matrix is given by 

j^2^( ml + hlHl{X,) hMtH2{X,)-^ilH,{X,))\ 
\ht{AlH2{X,) - fx,H,{X^)) ml + hlHl{X,) ) ' 

In thermal equilibrium the Green function for a quasiparticle with mass m is 
iA^{k^) = 27r5{k^-m^) sign{ko)fBE{ko) 

iA> (A;^) = 27rS{e - m') sign(A;o)(l + /se(A;o)) (6) 
with the Bose-Einstein distribution function 

The particle density can be deduced from the Green function using 

MX,) = 2zf ^k.A<{X„k,). (8) 

Since there will be already a contribution to the CP-violating particle densities in the 
mass term, we will in our approximation neglect interactions with other particle species. 
However, we will keep the collision term C, since this term usually drives the system back to 
equilibrium and allows to fulfill the physical boundary conditions far away from the wall. We 
will not explicitly calculate the collision term, but finally replace it by a phenomenological 
damping term. Hence the Kadanoff-Baym equations simplify to 



e 



-iO{k^ -M^}{A<'>} =C. (9) 



A further simplification is to perform the calculation in the bubble wall frame. Our picture 
of the phase transition is as follows. Bubbles of the Higgs field condensate nucleate and grow 
at a first order electroweak transition, and as they become large, they become approximately 
planar. The wall frame is then defined as the frame moving with the bubble phase interface. 
Due to the planarity, in this frame the mass matrix depends only on the average coordinate 
z :— Xs. In addition as mentioned in the introduction, we are working in the thick wall 
regime, what makes a gradient expansion reasonable. The system expanded up to first order 
in gradients reads (prime denotes derivatives with respect to z): 

i^k' + ikA + \dl -M'- \M^'dk}j A< = C. (10) 
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Using the hermiticity condition A^^ = — A"^ this equation can be spht into its hermitian 
and antihermitian parts 



k' + ^dl)A<-^{M^A<}- 







kAA< + [M', A<] - i {m", d,A< } = C 



(11) 
(12) 



where [■,■] and {■, ■} denote commutators and anticommutators. In the following we refer 
to these two equations as the constraint and kinetic equation. 



B. Lowest Order Solution 

Let us first discuss equations (fTTHT^ for a two-dimensional mass matrix that is constant 
in space and time. The mass matrix can be diagonalized by a unitary transformation and 
the equation in this basis reads (A^^ denotes the diagonalized mass matrix and A^ the 
corresponding Green function that is non-diagonal in general) 

(k' + \d!)^<-l{MlA<}=0 (13) 

kA^^ + '-[MlA<] =C,. (14) 

The question is, in which sense these equations can recover the solution in thermal equilib- 
rium (jni). We expect that the Kubo-Martin-Schwinger (KMS) equilibrium condition is then 
satisfied, such that Cd = 0. We can use the derivative of the second equation to obtain 

k'A< - [Ml [MlA<]] - \ {MlA<] = (15) 

kAA<+'-[MlA<]=Q. (16) 

Note that, upon the identification, iv)m {mm)) with A^^, these equations become identical 
to the leading^ order equations obtained for the chiral fermionic distribution functions 
(ql) in Ref. The constraint equation (fT3j) is algebraic, and it determines the spectrum 
of the quasiparticles in the plasma. At this point it is helpful to introduce two projection 
operators 

P^X = j-^[Ml[MlX]], P^ = l-P^, (17) 

where A := VTrAl^ — 4 DetAi'^ = Tr (a^Ail) denotes the difference of the eigenvalues of 
= I5 2, 3) are the Pauh matrices. The properties of the projection operators 
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can be easily checked. 

In the mass eigenbasis -P"^A^ corresponds to the complex off-diagonal entries, while 
P^A^ corresponds to the two real diagonal entries. If we split A^ in its transverse and 
diagonal parts A J := P^A"^, := A^ and using the relations 

{F^,X^} = 2r^X^, [Y^'.X^] = (Tr Y)X^ , (19) 

P^'Ml = Ml P^Ml = , (20) 
the constraint equations p5|) for the diagonal and transverse parts of A^ decouple 

(k' - Mi) Af = 0, (21) 

'''' 

Both diagonal and transverse constraint equation are algebraic, and thus the solutions are 
given by the appropriate 5-functions, which represent sharp on-shell projections. The di- 
agonal shell is given by the standard dispersion relation, whose frequencies are, fcg = uf = 
k"^ + mf, where mf are the eigenvalues of A^^. The transverse parts fulfill a different on- 
shell condition, which can be easily obtained from (j22j) . Note that these on-shell conditions 
are the same as the ones found in ^ by solving the leading order constraint equations for 
fermions. 

The kinetic equation (fTB|) reveals another difference between diagonal and transverse 
parts. The kinetic equations read 

hd^A^ = 0, (23) 
fcAAl + l[A<^,Aj] = 0. (24) 

The diagonal parts are constant in space and time, while the transverse parts rotate in 
flavour space with the frequency ~ A/kz- 

In the equilibrium solution (jUj) the transverse entries vanish everywhere, but it is clear 
that this oscillation dominates the dynamics of the transverse parts as soon as they are 
sourced by higher order contributions in the gradient expansion. 

Alternatively, oscillations can be induced by the initial conditions. This is, for exam- 
ple, the case in neutrino oscillations. Neutrinos are namely created as flavor eigenstates, 
and hence, from the point of view of the mass eigenbasis, a mixture of diagonal and trans- 
verse states. Since in most environments the damping of neutrinos is very small, neutrino 
oscillations persist for a long time. 
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C. First order solution and CP violation 

Let us consider again the Kadanoff-Baym equations (fTTHT^ to first order in gradients. In 
the last section we saw that in lowest order the spectrum can be separated into the diagonal 
and transverse contributions. One can show that in the first order system (fTTHT^ however, 
the different quasiparticles start to mix and the spectral functions acquire a finite width. 
This is reflected in the fact that, at flrst order in gradients, the constraint equation is not 
any more algebraic. 

Fortunately we do not need any information about the spectrum to solve the kinetic 
equation (fT2|) . since it does not explicitly contain any ko dependence. When transformed 
into the mass eigenbasis, the kinetic equation reads 

kAA< + [S, A<] + [Ml A<] - 1 {m'J + [S, Ml] , d,Al } = C, (25) 

with 

S = U^U' (26) 

and the matrix U{z) diagonalizes M'^, M"^ = WM'^U. 

The next step is to determine the CP violating contributions to the particle densities. 
By definition the CP conjugation acts as 

A^^(X, k) = CV A,{X, k) CV = A*(X, -k), (27) 
= (Xo, -X,), fc^ = (fco, -h) . (28) 

This transformation is in our equation (j25|) equivalent to 

U^U*, S^S*. (29) 

Now suppose that as in the chargino case our particles do not directly couple to the sphaleron 
process. Then the CP-violating particle density has to be communicated to the other species 
via interactions. Therefore we are rather interested in the CP violating densities in the 
diagonal matrix elements of the Green function in the interaction eigenbasis. These are 
given by 

Tr [A< - CVA<CV] = Tr [UA^U^ - U* Af^U^*] 

= Tr [U{A< - Af^nu'] (30) 
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and 

Tr [(73A< - (TgCP A<CP] = Tr [asUA^U^ - asU* Af^U^*] 

= Tr [asU iA< - A f-^nu'l (31) 

where the latter equahty in both cases follows from the fact that A^"^"^ is hermitian. Hence- 
forth we consider in the mass eigenbasis the equation for A^^ := A'^''^*. This Q-conjugation 
coincides with CP-conjugation on the diagonal, but it is in addition basis independent, since 
it commutes with the diagonalization matrix. This fact was already used in 1] to identify 
CP-violating quantities for mixing fermions before the Green function was transformed to 
the interaction eigenbasis. 

The equation for A^® is given by (notice that E is antihermitian) 

kAA<^ + K [S, A<2] - 1 [Ml A<S] - \[m'; + [S, Ml\ , d,At^] = C, . (32) 

The only change with respect to the original equation of A^ is a sign-change in the oscillation 
term [Ai^, A^*^] . If we include higher order terms in the gradient expansion additional Q 
breaking terms will appear. Since in leading order CP violation is based on the oscillation 
effect, one has to solve only the equation of the transverse parts and its Q conjugate. 
Collecting terms, that are at most first order in gradients (deviations from equilibrium 
6 Ad = All — Acq, M.'^' and S are counted as of order one in the gradient expansion) we get 
for the transverse deviations, 

kASA^+ '-[Ml6A^] -C, = S, 
KdJA^Q - [Ml 5 AJS] - = 5, , (33) 

with the source term 

S, = -h [S, A<]^ + i {mI + [J:,MI\ ,d,Af^Y . (34) 

This can be solved numerically using an Ansatz for a flow solution as described in Ref. 

Since A^'"'^ and A^*^ differ only by transposition, this calculation in addition shows that 
the diagonal entries in the mass eigenbasis will be CP even up to first order in gradients. 

III. THE DAMPING TERM 

If we solve equations ()33p without the collision term, we will have problems to ensure 
that our solution will be close to thermal equilibrium on both sides at a large distance from 
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the wall. This problem can be solved by introducing a damping term, that corresponds to 
statistical effects due to the interaction of the particles with the heat bath. In the rest frame 
of the plasma, the damping should take place in the positive time like direction as e.g. in 
the equation 

k ■ dxA + koTA = S (35) 

In the wall frame this leads to Cd — 7t)^(A;o — Vu,kz)rA with the wall velocity Vyj. For 
r a reasonable choice is F = aTc, where a denotes the couphng strength of the dominant 
interaction of the species and T^. is the temperature of the plasma during the phase transition. 

However, by introducing a term that breaks time reversal invariance, we run the risk of 
breaking CP explicitly by introducing new artificial CP-violating sources. We illustrate this 
by the following simple example. Assume that a quantity W, which denotes a CP-violating 
deviation from equilibrium, fulfills the equation 

d,W = exp(-z2) nBEiVW+r^) = ^i^^ K) (36) 

and we are interested in J dkz W{z,kz). 

To solve this equation, we can use the Green function method with the boundary condi- 
tion, such that W vanishes in the unbroken phase {z — > — oo), where the wall has not yet 
infiuenced the plasma. Then 

W{z, kz) = j dz' giz, z') S{z', kz) , (37) 

with the Green function 

g{z,z')^e{z-z') (38) 

and J dkW{z, kz) can be determined. 

Since the solution does not vanish in the broken phase {z +oo), we introduce a 
phenomenological damping term that breaks time invariance and our choice could be in 
analogy to the considerations above 

dzW + ^TW = S{z, kz) (39) 
kz 

The corresponding Green function is 

(ko/kz > : exp ( - {ko/kz)r{z - z')) e{z - z') 

9i^,y) = { / \ (40) 

[ko/kz < : -exp ^(A;o/A;^)F(^ - z')j 0{z' - z) 
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and yields the desired result. On the other hand, if the source is odd in k^, the picture 
changes. The equation (jHTj) yields a solution, that is odd in k^ and / dfe^ W{z, k^) vanishes, 
while the solution of gives a non- vanishing result even after integration over kz. 

The same effect can be seen in the kinetic equation Without the damping term, 

the result will be odd in kz-, such that only the three component of the particle current 

UX,) = 2il ^^k^A{X,,k,). (41) 

Jko>0 i^V 

is sourced. This is expected, since if this current is Lorentz boosted into the rest frame of 

,1 1 ji r^T~t • 1 J-' 1 J --L -plasma— frame -wall- frame • i 

the plasma, the ur-violatmg particle density Jq = 7^,^ Vyj vanishes m 

the static wall limit, = 0. 

After the damping term is introduced, j^°-^^°--f'^°-'^^ jg sourced even in the case of a static 
wall profile, which is clearly an unphysical result for a CP-violating quantity. Notice that 
this source persists even in the limit F/i — > 0. In the following we keep only the source terms, 
which are not induced by the damping term. 



IV. TRANSPORT EQUATIONS FOR MIXING CHARGINOS 

In this section we recall the fermionic transport equations derived in Due to the 
additional spinor structure of the Green function, we have to solve two equations for the 
left-handed and right-handed densities separately. In addition, the Green functions have 
a spin quantum number s. As in the bosonic case, only the transverse parts oscillate and 
contribute to the CP-violating (or better Q-violating) densities. In the mass eigenbasis the 
equations for 5g^^ and ^^f^^j read (see Eq. (78) in Ref. [l|) 

Kdz5g^^ + \ [ml + koTn5g^^ = S% (42) 
kASgl', + [ml 6gll\ + k,TJgll = SI , (43) 

with the spin-dependent part of the sources 

S% = -s^ [VV^\g,,,^] - -f [V{m^'m - m^m')V\ g,,,^] + 5^ {V{m^m)'V\go,e^Y 
ko 4A;o 4A;o 

SI = [f/f/f, ^o,eJ + [U{m'm^ - mm^')U\ go..^] - {U{mm))'U\ go,^^^. (44) 
ko 4:ko Ako 
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The function (^o.eq denotes the 70 coefficient of the Green function in thermal equilibrium 
and mass eigenbasis 

^o,eq = 2n\kQ\S{k'^ - mDfFD, Ifd = ^^^^^ ^ • (45) 



The chargino mass matrix m is given by 



m 




(46) 



and diagonalized by the biunitary transformation 

ma = UmV^ , (47) 

where /ic and M2 are the soft supersymmetry breaking parameters. 

To compare the result from these equations with the work 13 , 3| it is helpful to examine 



the contributions of the different sources in the local approximation, F/j ^ 00, in which 
diffusion transport is neglected. In this case, the resulting CP-violating vector and axial 
vector particle currents behave as 



Tr(a3j5,) =5,^ Tr(a3j,)=5; + 5^ 

= 2r-^53(M2/x,)(|M2p-|/ien9^(«i«2)r/fo) 
Si = 2T-'^iM2fic)iul-ul)d,iu,U2)vU 

S; = -2T-253(M2/ie) {u^d^u, - u,d^U2) (r/°o) + AT^f^^) , (48) 
where Ui^2 = 5'|-f^i,2|, and rj^Q^, t]^q^ and 77^2) are integrals derived in [l|. 



rr2-n I 1^ 1^0 , n \ 5 ""12/ c-/ 7 2 2 \ 



.o>o(2vr)3^o ' kiTl + {A/2y 
1 

= ^{VHl + V{n)2) , Vln) = ^{Vin)l ' V{n)2) , (49) 

and n = n{k^,m\2) denotes the distribution function. The contributions Sa and result 
from the ffist term in the sources ()44j) . while the term Sc results from the second and third 
terms in the sources (jSJ. 

Comparing with Eq. (3.13) of Q], we see that in local approximation our sources agree 
in the characteristics of the z-dependence, but show different structure in momentum space. 
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To facilitate a comparison with the work on semiclassical force baryogenesis of Chne, 
Joyce and Kainulainen we quote the dominating local source at the second order in 



re quot 

HQ. 



gradients in the plasma frame 

Tr (l jg) ^ = 2t;^T-^53(M2/ie)(w25^i + u,d',U2)Cl) ' (^0) 

where (^^^ = '?7^q)|a^o- This source corresponds to the CP-violating shift in the dispersion 
relation and dominates if the mass spectrum in the chargino sector is far from degeneracy 
(A —>■ oo) and in the limit of a small damping. It contributes in contrast to the first order 
terms to the trace of the chargino current. 



V. DIFFUSION EQUATIONS 



Using our formalism, we can deduce the CP-violating particle densities in the chargino 
sector. To evaluate the baryon asymmetry in the broken phase, we need to compute the 
density of left-handed quarks and leptons rip in front of the wall. These densities couple to 
the weak sphaleron and produce a net baryon number. 

To determine how the CP-violating currents are transported from the charginos to the 
eft-handed quarks and leptons we use a system of coupled diffusion equations as derived in 



3, 



and later adapted in 



13 



and 



v,„n 



10|. The diffusion equations are 
_r !^ - 

J- rn , 



kc 



v,„n 



w I I'M 



-er, 
+3r, 

Dhn'' + TY 



'nQ _ 






_kQ 


kp 


kn 






+ 9 


riQ + riT 




kp 


ks 




'riQ_ _ 


riT 


UH + rih 


_kQ 


kp 


kn 




riT 


+ 9 


Uq + Ut 




kx 


ks 




'riQ _ 


riT 


riH + rih 


_kQ 


kx 


ku 




'nQ_ _ 


riT 


riH + rih 


_kQ 


kx 


ku 





ko 



kp 



riT 
kp 



(51) 



J- h-r- 

kn 



:r. + 4r,,^ 



nh_ 

H 



(52) 
(53) 
(54) 



where ut denotes the density of the left-handed top and stop particles, uq the remaining 
left-handed quarks and squarks and nn and rih the sum and difference of the two Higgsino 
densities and nH2- The quantities fcj are statistical factors defined by rii = ki^i^ (/ij 
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denotes the chemical potential of species i). For light, weakly interacting particles ki ^ 2 
(bosons) or fcj ^ 1 (fermions), while for particles much heavier than T^, fcj is exponentially 
small. We use the values 

kQ^6, kT^9, kB^3, kH^l2 (55) 
corresponding to the light stop scenario j3l and the diffusion constants are j2l| 

~ 6/T„ ~ 110/T,. (56) 

B 



For the particle number changing rates we take 
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1111 
~ Yo'^'^' ^™ ^ 10 ~ 20'^'^' ~ 10'^"^ ^^'^'^ 



and for the sphaleron rates [2 

r,, ^ 1.5 X 10"2 T„ r^, ^ 6.0 X 10-^ T,. (58) 

The diffusion equations (jHTHH^ are derived under the assumptions that (a) the super- 
gauge interactions, which are of the weak strength, are in equilibrium; (b) the chargino 
asymmetry gets transported to the quark sector via the strong top Yukawa interactions, 
while the wino asymmetry does not contribute; (c) the gaugino helicity-flip interactions are 
in equilibrium, implying that the chemical potentials for particles and their supersymmetric 
partners are equal. These approximations imply that the main channel for baryon pro- 
duction is the conversion of the chargino asymmetry into the top sector, which then bias 
electroweak sphalerons. The accuracy of these approximations will be addressed elsewhere. 

The solution of Eqs. (j3TVI3^ is performed in several steps. First we use the transport 
equations in the chargino sector as described in [l| to determine tih and Uh- The result 
is used as an input in the equations (j51|) and (j52p . From these equations the left-handed 
particle density rii = dng + irix can be determined and used as a source for the weak 
sphaleron process as described in [3| (see also Ref. |3]). The net baryon density is given 
by 



^ ws 



and finally the baryon-to-entropy ratio is determined via 



/dz nL(z)exp \z ^ ^ \ (59) 



^ = ^, s = —g,^T^,^hl.lTl (60) 
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To check, whether our solution of the diffusion equation is consistent, we used the densities 
nq and rix as input for the equations (j^Hj) and (jMj) . The resulting deviations in the Higgsino 
densities never exceed 5% of the original densities. This is due to the fact that the Higgsino 
diffusion constant is rather large and that the oscillation partially suppresses an efficient 
transport of the quarks and squarks. In this light the equations of the Higgsinos decouple, 
since the oscillation provides the shortest time-scale. 

Note that in the work 2^ a suppression was found for the parameters of the Standard 
Model {kx ~ 3 in Eq. ()55p). As explained, for the mixing sources we consider here, the 
oscillation effectively decouples the dynamics of the charginos from the quarks/squarks. This 
allows us to neglect the backreactions from the quarks/squarks and leads to the absence 
of the suppression for kT ~ 3. If the oscillation is not the shortest time-scale, i.e. for 
l/ic — < 5 GeVjthe backreactions become large and our approach does not reproduce 
the suppression of 23| and would indeed over-estimate the result. In the following we will 
employ the parameters of Eq. (|55|) where this suppression mechanism is already ineffective. 



VI. NUMERICAL RESULTS 
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FIG. 1: The original Higgsino densities and the corresponding back-reactions. The parameters of 
the plot are iic = 200 GeV, M2 = 180 GeV, uia = 200 GeV 



In this section we will present numerical results of the transport and diffusion equations. 
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The Higgs vevs and the /3 angle are parametrized by 

H^iz) = H{z) sin(/5(z)), H^iz) = H{z) cos(/3(z)) (61) 

and 

H{z) = \v{T) [l - tanh (" " ) ) ' (62) 
m =P^- ^ A/3 (^1 + tanh (^a (^1 - ^ ^ . (63) 
The parameters used are 

Tc = 110 GeV, v{T) = 175 GeV, a = ^, tan{f3^) = 10, = 20/% (64) 
and the complex phase is chosen maximally 

3(M2/ie) = |M2/ie| • (65) 



We have checked, with the help of a program developed by the authors of Refs. [13|, |24| . 
that the values for v{T) compatible with present Higgs bounds typically lie in the range 
165-185 GeV. The exact value depends on parameters of the Higgs and squark sectors which 
affect our results only through this expectation value. We therefore have fixed the value of 
v{T) to its zero temperature result. The uncertainty arising from our choice is below ten 
percent. 

n 

The values of A/5 are deduced from [25] for the different values of rriA- The wall velocity 
is taken to be = 0.05 and the transport equations are evaluated using the fluid Ansatz 
for the first six momenta. The parameters of the diffusion equations are given in the last 
section. 

The plot Fig. supports the claim that, within our approximations and for our choice of 
parameters, the back-reaction of left-handed quarks and squarks, uq, ut, on the charginos 
can be neglected. The amplitude of the Higgsino densities coming from the back-reaction is 
always smaller than 3% and never leads to corrections of the baryon-to-entropy ratio larger 
than 5%. 

In Fig. 121 we plot the first order sources S^, S^, and the second order source (semi clas- 
sical force) 5°', as defined in Eqs. fl48fl3n|) . The first order sources are roughly of the same 
magnitude, and they peak when |/ic| ~ IM2I, where they also switch the sign. The second 
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FIG. 2: This plot shows the first and second order sources as a function of /ic with M2 = 200 GeV. 
The plot on the left are the sources with the damping, T = a^Tc, while on the right plot, T = 
{).2bayjTc. 

order source varies slowly with |/ic| and tends to dominates when the difference |/ic| — \M2\ 
becomes large. Note that when the damping is small, the first order sources become more 
peaked around |/ic| = IM2I, but the amplitude of the baryon asymmetry does not signifi- 
cantly change. On the other hand, the second order source is about an order of magnitude 
larger in the right plot, implying that in the limit of a small damping the second order 
source (semiclassical force) may result in a viable baryogenesis. Since our damping term 
is phenomenological and flavor blind, it would be premature to conclude that the second 
order source cannot lead to a viable baryogenesis, until a more quantitative analysis of the 
damping term is performed. 

The parameters chosen in Figs. El and 0] are similar to the ones chosen in jj^], in order to 
facilitate comparison. In plot Fig. El the parameters and M2 are fixed while /ic is varied. 
The maximum is not exactly at /ic = M2 as in but rather close to /ic ~ M2 + 20 GeV. 
The reason for this difference is that in our case all sources pH|l are of similar order, while 



m 



14| . the baryon asymmetry is completely dominated by a source term of form 5^ in 
that is proportional to A/5 in the parametrization (jUHj) and hence suppressed for large 
values of as shown in y. Another difference is that our plot shows the suppression 
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FIG. 3: This plot shows rjio = lO^^rj as a function of /Xc with M2 = 200 GeV and for several values 
of niA in GeV. 

for fic ^ M2 what is expected since in this case the quasi-particles have highly separated 
on-shell conditions and mixing should be suppressed. We would like to emphasize that the 
peak around /ic ~ M2 + 20 GeV is due to this suppression and not due to a resonance in the 
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sources as it was in the publications 
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26| and more recently in 
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In the present 



work, the sources show a resonance but the CP-violating densities do not since they are 
generated by the oscillations (see Eq. (j33|l ) and contain near the degeneracy an additional 
proportionality to the mass splitting A. In Fig. |3] the baryon asymmetry is plotted near the 
maximal value fic ~ M2 + 20 GeV. The maximum is reached near fic ~ 80 GeV in contrast 
to where the maximum was fic ~ 250 GeV. 

Finally in Fig. |S1 two contour plots are shown with regions in the (M2,/ic) parameter 
space for the baryon asymmetry expressed in terms of 7710 = lO^'' x 77. In these units the 
observed value is close to unity, 77io,obs = 0.8 — 0.9. If r^io > r^icobs, the observed value can 
be attained simply by adjusting the complex phase, which is in our calculation chosen to be 
maximal. The two plots correspond to the choices = 200 GeV and tua = 400 GeV. 

In the following we will comment on differences between the formalism used in this paper 
a„d the wo* H that lead to the d,serepa„ey ■„ the uumedea. results fl. As aheady 
mentioned in a previous section the authors in work in the flavour basis and write 
classical Boltzmann equations using CP-violating sources whereas in this work the sources 
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FIG. 4: This plot shows rjio = lO^^rj as a function of fic, M2 = /Lie — 20 GeV and for several values 
of niA (in GeV) . 

appear genuinely in a basis independent set of quantum transport equations. In this work 
the damping Th is primarily introduced to obtain consistent boundary conditions and it 
corresponds to the helicity flip rate Th in the diffusion equations of We have excluded 
Th dependent terms that violate CP and the limit F/i — > is straightforward. 

In addition to the damping Th a Breit-Wigner width ^^ was introduced in the chargino 
spectrum in We have checked in the simpler bosonic case that, for ^ 0, the present 



sources and those in jl3| are related in a simple way. A detailed discussion is presented in the 
Appendix. Of course, the ambiguity related to the magnitude of the source in the chargino 
diffusion equations remains in the formalism used in Refs. jisl-ll^. where a phenomenological 
thermalization time r or the classical Pick's law had to be used to incorporate the sources 
into the diffusion equations. In our formalism the magnitude of the source is completely 
specified. 

Furthermore, we have checked that the effect of the Breit-Wigner broadening on our 
sources is small. This effect can be modeled by replacing the ^-function in ()45p by the 
corresponding Breit-Wigner form. To account for the finite F^ in the transport and not just 
in the sources is on the same level as a treatment of the collision terms in Refs. 13, 
and it is outside the scope of this paper. In principle, the collision term could as well yield 
additional CP-violating sources, but a one-loop calculation [l^ in a model theory of chiral 
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FIG. 5: The baryon-to-entropy ratio tjiq = 10^^ x rj in the (M2,/Uc) parameter space from 
(0 GeV,0 GeV) to (500 GeV,500 GeV). For the left plot the value tha = 200 GeV is used, for 
the right plot rriA = 400 GeV. The black region denotes ??io > 1, where baryogenesis is viable. The 
other four regions are bordered by the values of tjiq, {—0.5,0,0.5,1}, beginning with the lightest 
color. 

fermions Yukawa-coupled to scalars, indicates that the collisional sources are phase space 
suppressed with respect to our tree level sources. 



VII. CONCLUSIONS AND DISCUSSION 



In this work we obtained the baryon asymmetry of the universe during the electroweak 
phase transition in the MSSM using semiclassical transport equations derived in a first 
principle approach from the Kadanoff-Baym (KB) equations in Ref. [3]. When the KB 
equations are expanded in gradients in the general case of mixing fermions, the CP- violating 
deviations from equilibrium can be sourced by a space-time dependent Higgs background 
both at first and second order in gradients. The first order effects, which occur only in 
the presence of fermion mixing, have been consistently determined including oscillations 
that are crucial for the dynamics of the CP-violating densities. The second order effects 
are dominated by the semiclassical force [7|, [101, 118|, which is the leading order source for 
single fermions coupled to a space-time dependent background. Unlike in some alternative 
approaches pursued in the literature, a nice feature of the present approach is that sources 
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and transport are treated within one formalism, which allows for an unambiguous fixing 
of the amplitude of CP-violating sources in (diffusion) transport equations. Moreover, this 
approach allows in principle for a systematic study of CP-violating sources from collisions, 
and how thermal and off-shell effects may affect the sourcing and transport of CP-violating 
charges. 

Furthermore, since our treatment is based on a formalism that fully includes the effects of 
mixing fermions, our results are manifestly basis independent. This is in contrast to former 
work, where the transport was treated either in mass eigenbasis [3, Q, Q|, or in flavour 
basis [3, ll^ ll^, and which describes just transport of two physical degrees of freedom, 
ignoring any dynamical effects arising from flavour mixing. For example, such a treatment 
of neutrino propagation would lead to complete neglect of neutrino oscillations. Unlike in 
the neutrinos case, the chargino oscillations occur on a microscopic scale given by the split in 
the chargino eigenvalues and by the chargino damping. In addition our formalism contains 
genuinely sources and transport such that no phenomenological thermalization time r has 
to be introduced as was done in [l3 . Il4 . Ilfil| . 

While a broad-brush picture of the first order sources resembles the sources found in 



Refs 
Nelson 



14| (this approach to supersymmetric baryogenesis was initiated by Huet and 



16|), there are noteworthy differences. Firstly, we found that chargino flavour oscil- 
lations are of crucial importance for identification and dynamics of the CP-violating sources. 
The oscillations tend to suppress the calculated baryon asymmetry, in particular in the limit 
of a moderate damping, a feature that was not observed in 14 1. Secondly, while broadly 
speaking the first order sources share similar parametric dependences with the earlier work, 
they do differ in some important aspects. 

Firstly, as can be seen in Fig. El all of the contributions to the BAU from our first 
order sources are of similar size, such that in the final BAU one sees the characteristics of 
all three sources. In particular, the BAU peaks at \fic\ — \M2\ + 20 GeV, and then dies 
out rather fast for large values of \fic\- On the other hand, the BAU obtained in Q| 
peaks at the chargino mass degeneracy, |/ic| — IM2I, it is about a factor 2 larger than in 
our calculation, and finally it does not diminish for large values of |/ic| as fast as in our 
calculation. Both discrepancies are due to the oscillations. Far from degeneracy (large mass 
splitting A) the fast oscillations will suppress the particle densities. Near the degeneracy 
(small mass splitting A) CP-violation is suppressed since it is generated by the oscillations 



22 



as shown in Eq. (jHHj) and this suppression cancels the resonance in the sources observed 
in Q. 

Provided it is not too strong, the phenomenological damping term that we introduce 
does not significantlv affect the maximum strength of the first order sources, unhke what 
wa. ob.e„ed . fllflQ. On the otKe. hand, the second o.de. soloes a. enhanced in 
the hmit of a small damping, as shown in Fig. |21 For example, for a moderate damping, 
r ~ a^T, the first order sources dominate in most of the parameter space. The second 
order source is small, such that that it cannot alone be a viable source for baryogenesis, 
even when the CP violation in the chargino sector is maximal. When damping is weak, 
r ~ 0.25awT, the second order source dominates in a large section of parameter space. For 
even smaller values of F the semiclassical force source alone represents a viable baryogenesis 



source, implying that our source is somewhat larger than what was found in Ref. |10|, which 
agrees quite well with the BAU found in based on a study of semiclassical force source 
obtained in the mass eigenbasis [isj ]. 

Perhaps the most severe constraints on the supersymmetric baryogenesis in near future 
are expected from electric dipole moment (EDM) measurements. Already the current con- 



place rather strict constraints on the CP-violat 



mg 



straints on the EDM of the electron 

phases in the chargino mass matrix, as can be seen from Fig. 4 in Ref. 0] or Fig. 6 in 
that claims a little less restrictive bounds. For example, for fic = 200 GeV, Mh+ = 170 GeV 
and tan(/3) = 6 the CP-phase is restricted to be less than about 1/12 and 1/10, respectively, 
implying that, when our numbers are taken at the face value, the baryogenesis mechanism 
presented here is by about factor 5-6 too weak to account for the observed BAU. Similar 
conclusion is reached for other values of |/ic| and IM2I since both the EDMs and the pro- 
duced baryogenesis decreases with decreasing chargino masses. We would like to emphasize 
that most of the parameters are chosen in order to produce as much baryon asymmetry as 
possible, e.g. the values used for the wall velocity Vw and the wall width Lyj. The only 
relevant parameter we have not varied so far is tan/5. Smaller values of tan/3 lead to less 
restrictive EDMs and at the same time to more baryon asymmetry as shown in Fig. IHl On 
the other hand tan (3 due to the mass of the lightest Higgs is restricted to be in the range 

5 < tan /5 (66) 
such that we are not allowed to enter the region with smaller values of tan/3. In addition 
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FIG. 6: This plot shows rjiQ = W^^rj as a function of /ic, M2 = /ic — 20 GeV, tua = 150 GeV and 
for several values of tan/3. 

even for tan/3 = 3 our result is still a factor 2 to small. Note that this is a very different 
conclusion from the one reached in Ref . , where an ample region of parameter space was 
claimed to result in a successful baryogenesis in the MSSM. 

Based on our analysis, can we conclude that the MSSM baryogenesis is ruled out? At 
least a factor 2 can be accounted for based on the inaccuracies in the parameters in diffu- 
sion equations, as well as from approximations that lead to the set of diffusion equations 
considered here, but unlikely a factor 5 Nevertheless, it would be premature to claim 
that the MSSM baryogenesis is ruled out, since the chargino mediated baryogenesis studied 
here does not exhaust the possibilities of the MSSM. Recall that neutralinos mediate baryo- 
genesis as well, and that their contribution may be as important or even more important 
than that from charginos. Furthermore, in the complete set of diffusion equations, there 
may be additional channels, which lead to baryon production, as of yet unaccounted for. 
Finally the EDM analysis given here is not conclusive. For larger values of rriA and due to 
the possibility of fortuitous cancellations between different EDM contributions, the value of 

n 

the electron EDM could be reduced relative to the generic values used above |2S|- Hence 



electroweak baryogenesis seems to 



with the conclusion drawn in 



m 



oe still possible in the MSSM and in this respect we agree 
14, ^|. However, we would like to emphasis two robust 



and novel consequences resulting from the quantum treatment of transport in the chargino 
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sector: The BAU is strongly suppressed away from the chargino mass degeneracy and one 
requires even close to the degeneracy a CP-violating phase of order unity, more precisely 
arg(M2/ie) > 0.2. 

Modifications of the MSSM with an additional singlet (NMSSM) also contain the promi- 
nent chargino channels. As shown in j^^, 2^ one easily can get a strong first order phase 
transition and also spontaneous CP-violation at the temperature of the phase transition 
not affecting the EDM at T = 0. This then allows for a satisfying baryon asymmetry with- 
out squeezing the (unfortunately many) parameters. One can also think about extensions of 



the MSSM that not forbid tan/? ~ 1 or modifications of the Standard Model where 
the chargino system does not appear, but of course again quantum-transport is important. 

In summary, our numerical solution to the diffusion equations shows that a 

successful baryogenesis at the electroweak scale from charginos of the MSSM is possible only 
when CP violation is quite large, and near the resonance, |/ic — M2I <^ 50 GeV, M2, ^ 
500 GeV. As long as the first order sources dominate, due to the oscillations, this picture 
persists also for much stronger sources, which is to be contrasted to 14. 15^. 

Our conclusion is that, in purely chargino mediated MSSM baryogenesis the capability 
to explain the observed baryon asymmetry is strongly constrained by the current electron 
EDM bounds. 
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APPENDIX A: COMPARISON OF BOSONIC SOURCES 



In this Appendix we show how the sources, presented in the current work, relate with 

nn □ 

those of references |13L I14| in the bosonic case and in the limit of zero widths . In order to 
inspect this we make use of the Kadanoff-Baym equations for the full 2x2 Green functions 
of the Schwinger-Keldysh formalism. These equations are obtained from Eqs. (jllllT^ if we 
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substitute A"^ by the full 2x2 propagator, 

insert unity in the r.h.s. of and set the collision term of ()12|) to zero, 

{k' - \d'x) A - I {M^ A} + 1 A] = 1 (A2) 

k-dxA+'- [M\ A] + i {d^M\ d^,A} = . (A3) 

In the corresponding Dyson-Schwinger equation for A is iteratively solved as an 
expansion in powers of d^Ai'^, 

A = A(°) + A(i) + . . . (A4) 

where A^^^ is the leading order equilibrium propagator, and A^^^^ = 0{d^A4'^) denotes 
a first order correction. Upon performing a Wigner transform over the spatial variables, 
J d^{x — y) e*(^~^)''^, identifying z = (x + ?/)/2 = X, and transforming into the flavour basis, 
the first order correction, A*^^) = A*^^''(A;;X) given in Eq. (2.5) of Ref. becomes 

A« = - [(a,^A(°))(ajA^2)A(°) - AW(9JA^2)^m^(o)] _ (A5) 

,n ..e app.oacK a„ . Q . ca,cu«on of .... one . — .n 

long range effects, and hence the term d\/S.^^^ in the constraint equation ()A2|1 and k ■ dx^'^^^ 
in ()A3|) were considered of second order, and thus have been neglected. From Ref. ill] and 
this work we know however that, when the dynamics is taken account of, in the case of 
mixing scalars and fermions the flavour oscillations mess up the derivative expansion, such 
that only the terms containing spatial derivatives acting on the mass term are genuinely 
der i vat i ve- suppressed . 

a on inserting ()A4|) into ()A2HA3|) and using the prescription for derivative counting 
, we get for the leading order propagator, 

A;2AW-i{A^2^AW} = l, (A6) 

which is solved by the thermal Green function, which commutes with M.^. The first order 
equations are, 

fc2^(i) _ i [M\ A(i) } + [d^M\ 9,^AW] = (A7) 
A;■9xAW + ^[^^A«]+^{9^^<^9,^AW} = 0. (A8) 
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By a judicious use of ()A6|) and its derivatives, 

- A^2)a^A(°) = -(2A:^)A(°) = (9^A(°))(/t2 - A^') (A9) 
one finds that the first order correction (IA5I) can be recast as, 



A(i) = ^[(9,^A(°))(A;2-A^2)^x^(o)_(^x^(o))(^2__^2)^MA(o)] 

= -i [A(°)A; ■ 9xA(°) - {k ■ 9xA(°))A(°)] . (AlO) 

It can be easily shown that, when this is inserted into ()A7HA8|) . one obtains two consistent 
equations for A^^\ 

Note that taking moments of the kinetic equation (IA8I) a llows for a simple prescription 
on how the CP-violating source originally calculated in [131] enters the relevant transport 
equations for squarks. The term A*^^) enters through the commutator [tVI^, A^^^] in ()A8jl . 
while used a heuristic prescription for the sources based on the Pick's law and interpreted 
the diagonal entries of A^^-* in the interaction basis as sources for the classical diffusion 
equations. 



Note further that, even though we have rephrased the source of in our language, it 
remains a nontrivial matter to establish the exact correspondence between the source of 
appearing in ()A8j) and the source calculated in this work. Our source is in principle obtained 
by the means of the kernel of Eq. acting upon (jBlI), which is thus of a complicated non- 
local form, and bares no simple relation to the source in ()A8|) . apart from a rather superficial 
similarity, stemming from the fact that the kernel of Eq. (jHHj) is a nonlocal functional of the 
commutator [tVI^ (/),-] acting upon (j3i|) (see Ref. Q). 

Pinally, we emphasize that the difference in how the sources couple to the diffusion 
equations cannot alone explain a different baryon asymmetry obtained by the two methods, 
but also the presence of the oscillatory terms. 

As regards the case of mixing fermions, we expect that the sources can be related in a 
similar fashion. Because of the spinor structure however, the comparison for fermions is a 
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nontrivial generalization of the bosonic case. 
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